Mon. Not. R. Astron. Soc. 000. ITHT81 (20081 Printed 24 August 2009 (MN M£X style file v2.2) 



O 
O 

(N 

6 

o3 



> 

OO 

cn 
o 

On 
O 



13 



Precision Astrometry of a Sample of Speckle Binaries and 
Multiples with the Adaptive Optics Facilities at the Hale 
and Keck II Telescopes 

K. G. Helminiak 1 *, M. Konacki 1 ' 2 , S. R. Kulkarni 3 and J. Eisner 4 

1 Nicolaus Copernicus Astronomical Center, Department of Astrophysics, ul. Rabianska 8 , 87-100 Toruri, Poland 
2 Astronomical Observatory, A. Mickiewicz University, Sloneczna 36, 60-286 Poznan, Poland 

3 Division of Physics, Mathematics and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 
4 Steward Observatory, University of Arizona, Tucson, AZ 85121, USA 



Accepted . Received ; in original form 



ABSTRACT 

Using the adaptive optics facilities at the 200-in Hale and 10-m Keck II, we observed 
in the near infrared a sample of 12 binary and multiple stars and one open cluster. 
We used the near diffraction limited images of these systems to measure the relative 
separations and position angles between their components. In this paper, we inves- 
tigate and correct for the influence of the differential chromatic refraction and chip 
distortions on our relative astrometric measurements. Over one night, we achieve an 
astrometric precision typically well below 1 miliarcsecond and occasionally as small as 
40 microarcseconds. Such a precision is in principle sufficient to astrometrically detect 
planetary mass objects around the components of nearby binary and multiple stars. 
Since we have not had sufficiently large data sets for the observed sample of stars to 
detect planets, we provide the limits to planetary mass objects based on the obtained 
astrometric precision. 

Key words: astrometry - planetary systems - instrumentation: adaptive optics - 
binaries: visual - stars: individual(56 Per, GJ 195, GJ 300, GJ 352, GJ 458, GJ 507, 
GJ 569, GJ 661, GJ 767, GJ 860, GJ 873, GJ 9071). 



1 INTRODUCTION 

In the field of exoplanets, astrometry has not had too many 
triumphs to date. Of over 350 known planets or plane- 
tary candidates only one has been discovered astrometri- 
cally jPravdo fc Shaklanll2009l ). The astrometric results pre- 
sented by Han et alj 1 200 J ) are disputable and the true 



masses of only a few planets were calculated more reliably 
by combining the radial velocities (RV) and astrometry from 
the H ubble Space Telescop e and the ground-based observa- 
tions (Benedict et al. 2003. [20061 ) . Nevertheless, astrometry 
may just turn out to be the most promising planet de- 
tection me thod in the future . Astrometr ic space m i ssions , 
like SIM (|Unwin et all 120081 ) or Gaia (jPerrvmanl 120051) . 
and a fe w ground-based interferom etric surveys are ongo- 



ing (e.g. lLane fc Muterspaug 



on the VLTI: 



Eisenhauer et al 



2004 ) or in preparation (e.g . 



20081 : lLaunhardt et al.ll200Sl ; 



ISahlman et aL 20081 ). In particular ground based interfer- 



ometers seem to be already well suited to detect planets by 
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providing microarcsecond (pas) astrometric precision (e.g. 
iMuterspaugh et al.ll2005T ) for bright nearby binary stars. 

The milliarcsecond (mas) or better precision can be 
achieved by imaging with the adaptive optics (AO) systems. 
This was already demonstrated for two bi naries, HP 19063 
and HP 19994 obser ved with the VLT (|Neuhauser et all 
l2006l : iRoll et alj|2008t ) a nd a globular cluster M5 observed 
with the Hale telescope l|Cameron et alj|2009l ). Such a pre- 
cision can be reached by means of relative a strometry over 
a small field of view fcameron et al l I2009D . To this end 
one needs to have at least one reference object not too far 
from a science object. For this reason visual and speckle 
binaries become a natural target for such measurements. 
Incidentally, the subject of the existence of exoplanets in 
binary and multiple stars has become of significant in- 
terest (e.g. Raghavan et al.ll2006l: lEggenberger et al.l [20071 : 
IMuterspaugh et all 120071 : iMugrauer fc Neuh &uscr 20091) . It 
is now accepted that the detection or lack of planets in star 
systems will provide add itional constraints to our models 
of planet formation (e.g. IHolman fc Weigert] 1 19991 : INelsonI 
l200d : ILissauer et alll2004[ ) . 
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Table 1. Number of images of a particular object per night. NGC 68*7 1 



Night /Tel. 


56 Por 


GJ 195 


GJ 300 


GJ 352 


GJ 458 


04 Mar./K 


58 


— 


58 


— 


— 


23 Apr./H 


— 


— 


— 


53 


975 


9*1 Tun 

j un . / n 










1060 


24 Jun./H 


— 


— 


— 


— 


685 


11 A >> t /TU 

zl Aug./n 




:;i h } 








22 Aug./H 




582 








13 Nov./H 





049 











Night/Tel. 


GJ 507 


GJ 569 


GJ 661 


GJ 767 


GJ 860 


04 Mar./K 




29 








23 Apr./H 


949 




656 






23 Jun./H 


1012 




454 




189 


24 Jun./H 


520 




800 




1166 


26 Jun./H 






1250 






21 Aug./H 






750 


509 


000 


22 Aug./H 






636 


746 


507 


13 Nov./H 








745 


584 


Night/Tcl. 


GJ 873 


GJ 873B 


GJ 9071 




NGC 6871 


23 Jun./H 


225 


251 






510 


24 Jun./H 


510 


497 






1010 


21 Aug./H 


20(1 


200 


750 




1083 


22 Aug./H 


200 


200 


513 




2131 


13 Nov./H 


300 




1246 




624 



A faint, third component of the GJ 860 system was not always in the 
field of view (due to the dithering). Separate columns for GJ 873 and GJ 
873B arc to distinguish between the images of the full triple system and 
the double secondary only. 'K' stands for Keck II and 'H' for the Hale 
telescope. 



In this paper, we present our observations of a sample 
of 12 binary and multiple stars and one open cluster ob- 
tained in 2002 with the AO facilities at the Hale and Keck 
II telescopes over the period of 7 months. We investigate the 
influence of several systematic effects that have an impact 
on the relative astrometry and demonstrate that by correct- 
ing for them one can achieve a sub mas precision. Finally, 
we provide the limits to planetary mass objects around com- 
ponents of our target stars derived from the obtained astro- 
metric precision. 




Figure 1. A field in NGC 6871 centreed on a = 20 ,i 05 m 57 s , 8 = 
35°47'25" as seen on 23 June 2002. Positions of only 5 labeled 
stars were measured. Other 3 (top of the image) are not seen in 
images from other nights. FOV is 25x25 arcseconds. North is up, 
East is left. 



used a 1%-transmission neutral density filter ND-1 for de- 
creasing a flux from very bright stars. Dithering was carried 
out by shifting the observed position by ~ 2 as. 

We also had one clear night at the 10-m Keck II tele- 
scope. Using its AO system and NIRC2 (the Near InfraRed 
Camera 2) we obtained data for three targets and the total 
of about 600 images. NIRC2 is a mosaic of four 512 x 512 
InSb Aladdin-3 detectors. For the observations we used 9.942 
and 39.686 mas/pix scale and the J, K' and K-cont (narrow- 
band) filters. Dithering was done using field rotation. 



2 OBSERVATIONS 
2.1 Instrumentation 

The main instrument in our project was the 200-in 
Hale Telescope at the Palomar Observatory. We used 
PHARO (the Palom ar High Angular Resolution Observer, 
lHavward et al.l lioOll ) camera with PALAO (the PALomar 
Adaptive Optics) system. PHARO uses a mosaic of four 
512 x 512 HgCdTe H AWAII detectors for observations be- 
tween 1 and 2.5 /im l|Havward et al.ll200li ). PALAO is an 
adaptive optics system mounted at the Cassegrain focus of 
the telescope. It employs a Shack-Hartman wavefront sen- 
sor and a Xinetics Inc. 349/241 active-element deformable 
mirror. PALAO 's detailed description can be found on- 

lineE 

With Hale/PALAO we obtained about 30,000 images of 
our targets. The data were collected over 7 nights nights be- 
tween April and November 2002. We used an imaging mode 
with 39.91 and mostly 25.10 mas/pix scale and K, K', Ks 
broad, as well as Br7 and Fell narrow band filters. We also 



2.2 Objects 

From all the attempted objects we selected 9 bina- 
ries/multiples from the Hale data set and 3 multiples from 
the Keck II one. These are: GJ 195, GJ 352, GJ 458, GJ 
507, GJ 661, GJ 767, GJ 860, GJ 873 and GJ 9071 for the 
Hale, and GJ 3OC0, GJ 569 and 56 Per for the Keck II. All 
systems are shown in Figure [2] We also selected a field in the 
open cluster NGC 6871 (the Hale data set) centreed around 
a = 20 h 05 m 57 s ,<5 = 35° 47' 25" as a reference field to study 
the systematic effects (Figure [T}. 

The selection criterion we adopted for the Hale sample 
was mainly the high number of individual images and also 
the number of nights during which a given object was ob- 
served. The exception is GJ 352. It was selected to study a 
precision of astrometry for a low number of single images. 
For the Keck II observations it was mostly important to 
check how good the AO correction was and how many un- 
saturated stars were in an image (see Figure [2}. Due to sat- 
uration we were able to measure the relative positions only 



1 http:/ /ao.jpl. nasa.gov/Palao/PalaoIndex. html 



2 Two fainter stars seen close to GJ 300 and GJ 873 are actually 
field stars. 
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Figure 2. Images for all observed systems and labels of investigated stars. North is up and East is left except for 56 Per where North 
is up-left (rotation by 44.7 degrees, as shown). FOV size varies, separations are given in Table [6] The saturated cores of the primary 
components of 56 Per and GJ 300 are clearly seen. GJ 569 A is also saturated but this is not seen in the adopted scale. 
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Figure 2 - continued 
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Table 2. Basic information on our targets. 



Star 


No. 


Sp. Type 


Magn. 


(Band) 


7r [mas] 


Comment" 


Telescope 


Ref. 


56 Per B 


1+2 


??? 


8.7 


(V) 


24.00(.91) 


double 


Keck II 


1.2 


GJ 195 A 


1 


Ml 


10.16 


(V) 


72.0(.4) 




Hale 


3,4 


GJ 195 B 


2 


M5 


13.7 


(V) 


72.0(.4) 




Hale 


3,4 


AG+45 517 


3 


77? 


11 


(V) 


??? 


field 


Hale 


4 


GJ 300 B 


1+2 


K7III? 


8.39 


(•1) 


125.60(.97) b 


double, field 


Keck II 


5.6 


GJ 352 A 


1 


M4 


10.07 


(V) 


94.95(4.31) 


— 


Hale 


1.7 


GJ 352 B 


2 


M4 


10.08 


(V) 


94.95(4.31) 




Hale 


1.7 


GJ 458 A 


1 


M0 


9.86 


(V) 


65.29(1.47) 


— 


Hale 


1.8 


GJ 458 B 


2 


M3 


13.33 


(V) 


65.29(1.47) 


— 


Hale 


1.8 


GJ 507 A 


1 


MO. 5 


9.52 


(V) 


75.96(3.31) 


— 


Hale 


1 


GJ 507 B 


2 


M3 


12.09 


(V) 


75.96(3.31) 




Hale 


1 


GJ 569 Ba 


1 


M8.5V 


11.14 


(J) 


101.91(1.67) 




Keck II 


1,9,10,11 


GJ 569 Bb 


2 


M9V 


11.65 


(J) 


101.91(1.67) 




Keck II 


1,9,11 


GJ 661 A 


1 


M3 


10.0 


(V) 


158.17(3.26) 




Hale 


1,7 


GJ 661 B 


2 


M4 


10.3 


(V) 


158.17(3.26) 




Hale 


1.7 


GJ 767 A 


1 


Ml 


10.28 


(V) 


74.90(2.93) 




Hale 


1.8 


GJ 767 B 


2 


M2 


11.10 


(V) 


74.90(2.93) 




Hale 


1,8 


GJ 860 A 


1 


M3 


9.59 


(V) 


249.53(3.03) 




Hale 


1,12 


GJ 860 B 


2 


M4 


10.30 


(V) 


249.53(3.03) 


flare 


Hale 


1,12 


CCDM 22281. ..H d 


3 


??? 


13.8 


(V) 


??? 


field 


Hale 


13 


GJ 873 A 


1 


M3.5e 


10.09 


(V) 


198.07(2.05) 


flare 


Hale 


1 


GJ 873 B 


2+3 


G 


10.66 


(V) 


198.07(2.05) 


double, field 


Hale 


1,14 


GJ 9071 A 


1 


K7 


10.2 


(V) 


72(4) 




Hale 


1,8,13 


GJ 9071 B 


2 


MO 


14 


(B) 


72(4) 




Hale 


1,13 



a If 'double', magnitude refers to a total magnitude of both components and spectral type is 'averaged'. If 'field', the star is not 

gravitationally tied with brighter components. 

b Parallax is for GJ 300, not the two investigated stars. 

c ISimon et ahl l|2006h suggested that GJ 569 Ba may be a binary with similar brightness components. 

d CCDM 22281. ..H = CCDM J222 81+5741H - a part of a m ulti-stellar system wh ich includes also GJ 860. 

Ref.:(l) The Hipparcos C atal ogue dPerrvman et al.lll997h: (21 iBarstow et al.ll2005l: C31 Ijenkindll952l : (4) The PPM No rth Catalogue 
iRoeser fc Bastian|[l988h; ( 5l lSimons et al.lll996t ( 6l" lrlenrv et al.ll2006l; (7 1 lAl-Shukri et al.lll996l ; (8l lReid et~alll2004t (9) 
lLane fe MuterspaugbfliixM (lOllSimon et al.ll2006l: (l ll lCutri et al.ll2003l: (121 lLaw et al.ll2008t (131 CCDM - Catalog of Components 
of Double & Multiple stars JPommanget fc Nvdl2002l) ; ( 141 lOppenheimer et al.ll200ll 



for the double secondary components of the Keck targets. 
The last two issues were caused by the varying weather con- 
ditions and hence highly variable AO correction. The final 
numbers of individual images per night used in the analysis 
for a particular object, after rejection of useless data, are 
given in Table [T] 

Our objects are mostly M-type dwarfs located less than 
20 pc from the Sun. In a few cases (e.g. GJ 195) not only 
a binary but also other stars were captured in an image. 
Their relative positions were also measured. For the open 
cluster, the astrometry was performed only for 5 stars but 
other objects can be seen in images as well. For all the stars 
which are cataloged the basic information is given in Table 
[2] Column 'No.' refers to a number/label in an image. The 
higher is the number assigned, the fainter the star is. 



3 POSITION CALCULATIONS 

The images were first reduced with standard IRAF tasks 
for data reduction. Subsequently, given the number of expo- 
sures, the relative positions of the stars were computed with 
our own software in an automated manner as follows: (1) 
the shifts from image to image (due to the dithering) were 
measured by cross-correlating the template image (usually 



the first image) with all the subsequent exposures and the 
approximate positions of stars in a given image were calcu- 
lated with an accuracy of ±3 pix (~ 75 mas for most of the 
Hale's sample) , (2) based upon these positions, the centroids 
were calculated, (3) based upon the centroids, the following 
2-dimensional elliptical gaussian function was fitted to the 
cores of the images of stars: 



G(x, y) = B + Aexp 



[(x — xg) cos 0-{y-yo) sin t 

[(x-xg) sin B+(y-y ) cos S} 2 
2^ 



(1) 



where B is a background level, A the amplitude of a gaus- 
sian, (xo,yo) is a position of the star, a x ,o- y the corre- 
spondimg widths and 9 is a tilt of the gaussian. 

We have decided to use such an approach because it 
offers a simple and robust way of modeling of the cores of 
stars' images. One could envision using an empirical point 
spread function (PSF) as a model for the images of stars. 
This is however challenging due to a fact that in a single 
image we typically have only two stars. Hence our knowl- 
egde about the actual empirical PSF for a given frame is 
limited. Additionally, since the PSF's shape varies, it is not 
practical to use several subsequent exposures as a reference 
for an averaged empirical. This is demonstrated in Figure 
[3] for a series of 5 frames of GJ 661 taken on June 23 and 
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Figure 3. Top left: A combined image of GJ 661 from five exposures taken on June 23. The corresponding frame numbers are next 
to each image of the binary; Top right: A residual image after a subtraction of an empirical PSF (computed with DAOPHOT; Stetson 
1987). The size of the fitting aperture (9 pixels in radius) is shown as a white circle. Upper middle: A 15 X 15 pixel zooms on the contour 
plots of the fitted gaussian functions at the positions of respective stars. Contour levels vary. Gaussians clearly show significant variations 
of their shape; Lower middle: The same zooms on the residuals after a subtraction of the gaussian functions. Changes of the first Airy 
ring can be seen. Color scale is the same in every sub-panel. Bottom: The same zooms on the residuals after a subtraction of an empirical 
PSF. A clear leftover in the PSF core is evident. Color scale is the same in every sub-panel. 



spanning 15.5 minutes for which an average empirical PSF 
is calculated using a 9 pixel aperture and subtracted from 
the images of stars. As can be seen a fitting of the gaussian 
performs better. The details of this procedure are also given 
in Table El 

The results of the gaussian fitting were used to compute 
the relative separations and the position angles of pairs of 
stars. The NIRC2 data were corrected for the field rota- 
tion used for dithering. Let us note that we did not use any 
weighting scheme for individual images as it was done by 
ICameron et al.l (|2009h . A significant improvement in the as- 
trometric precision after using an optimal weighting is seen 
when the number of reference stars exceeds 5 (see Fig. 2 in 



ICameron et aT]l2009T) which is never the case in our work 
targeting binary stars and aimed at investigating the astro- 
metric precision in the case of only one reference star. 



4 SYSTEMATIC EFFECTS 

4.1 Adaptive Optics correction and field of view 

The main factor allowing us to obtain precise astrometric 
measurements is obviously the adaptive optics. Its perfor- 
mance will influence the final astrometric precision as can 
be seen in the case of GJ 352 (p ~ 350 mas) observed dur- 
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Table 3. Parameters of an elliptical gaussian function used to model the images of stars in the frames for GJ 661 (Figure [3}. An average 
value of a residual after a gaussian and an empirical PSF subtraction is given together with its uncertainties in the last four columns. 
Analyzed frames are numbered as in Fig. \3\ 



Star 


A 


(fx 


<Ty 


x 


yo 


B 


e 


Av. resid. 


(j 


Av. resid. 


<7 


Frame/ID 


[counts] 


[pix] 


[pix] 


[pix] 


[pix] 


[counts] 


H 


(gauss) 




(PSF) 




No. 50/1 


2846.443 


1.505 


1.432 


543.897 


615.032 


13.231 


130 


17.59 


18.16 


-111.73 


13.57 


No. 50/2 


2203.045 


1.437 


1.429 


548.635 


588.264 


48.128 


91 


3.74 


11.38 


-101.62 


10.73 


No. 100/1 


2711.196 


1.399 


1.499 


622.879 


694.661 


-0.042 


122 


4.88 


11.46 


-39.33 


5.68 


No. 100/2 


2125.978 


1.378 


1.512 


627.583 


667.851 


-0.300 


125 


-3.95 


14.32 


-85.35 


7.45 


No. 150/1 


2629.679 


1.478 


1.407 


626.141 


538.949 


35.679 


16 


11.13 


11.38 


-130.78 


10.17 


No. 150/2 


1953.843 


1.636 


1.527 


630.955 


512.050 


-17.370 


65 


15.58 


17.12 


-116.19 


17.75 


No. 200/1 


2391.502 


1.395 


1.488 


462.866 


535.715 


59.381 


115 


3.06 


10.35 


-140.15 


8.43 


No. 200/2 


1928.833 


1.363 


1.469 


467.544 


508.818 


51.934 


65 


8.88 


13.26 


-75.02 


7.67 


No. 250/1 


2344.759 


1.436 


1.349 


626.177 


652.063 


81.791 


165 


2.47 


7.87 


-112.96 


9.03 


No. 250/2 


1722.767 


1.441 


1.450 


630.861 


625.335 


50.023 


87 


0.83 


6.09 


-67.74 


7.36 




100 120 140 

Image number 




100 120 140 

Image number 



Figure 4. Impact of Adaptive Optics correction on relative position measurements in case of GJ 352. After image No. 107 (dashed line) 
AO works improperly. 



ing challenging weather conditions. From all 75 images of 
GJ 352 taken only the first 10 were properly corrected (the 
Airy pattern visible) and only in 53 images components were 
resolved and could be analysed. For these 53 images the 
centroids were calculated. Subsequent gaussian fitting was 
possible only for 34 images for which the fitting procedure 
converged. The centroids and the outcome of gaussian fitting 
are in agreement for the first 10 (Figure |3J. 

Another factor having an impact on the astrometric pre- 
cision is the field of view and the corresponding PSF sam- 
pling which was especially important for the Keck II data. 
The images were taken in two pixel scales (field sizes) - 
9.942 mas/pix (10" x 10" field, narrow) and 39.686 mas/pix 
(40" x 40" field, wide) and in various field rotator angles - 
45.7, 0.7 and —44.3 deg. For this telescope the diffraction- 
limited size of a star's image in the near infrared corresponds 
to about 1.4 pixel in the wide field. For faint stars it means 
that most of their light is collected in one pixel which makes 
the PSF undersampled and the gaussian fitting difficult. 
This issue, however, can be at least partially overcame by a 
sub-pixel dithering. For the purpose of this paper, we used 
only the frames taken with the 'narrow' field. 



4.2 Atmospheric refraction 

The atmospheric differential refraction (ADR) creates a shift 
of star's image. It is highly dependent on the zenith angle 
and wavelength. Formulae for computing ADR effect are 
given for example by Roe (2002) where the angle R, which 
is the difference between the real and observed zenithal dis- 
tance is given by: 

R = z t - z a ~ 206265 0^2~ j tan Zt [ arcsec \ ( 2 ) 

where z t is the true zenithal distance, z a is the observed 
zenithal distance, and n is the refraction index, dependent 
on the wavelength A and weather conditions: 

n{\,p,T,p w ) = 1 

+ [64.328 + T |e^ + 3 a^]^10- ( 3 ) 
_43.49 [1-24°wm] at 10 -« 

where A is given in fim, p, T and p w are the pressure [hPa] , 
temperature [K] and partial pressure of water vapor respec- 
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tively. Symbols with the index s refer to the canonical values 
of air pressure (1013.25 hPa) and temperature (288.15 K). 
The angle R is much smaller in IR than in visible. 

ADR also affects relative astrometric measurements. 
Since two objects are seen at different zenithal distances 
Zi and Z2, the corrections Ri and R2 are also different. 
The component of the separation vector parallel to the 
direction to the zenith increases after ADR correction by 
AT? = [iZ_ — Ri\. This quantity changes with weather con- 
ditions ( air pressure and temperature). As we have demon- 
strated (|HeIminiak| [2009), the magnitude of this change is 
often higher than an achievable astrometric precision even 
for relatively compact systems. Clearly, ADR's influence 
must be corrected for and the weather conditions should 
be well known. This c onclusion is in contrary to that by 
iNeuhauser et all (|2006T ). who claim that the refraction is in 
general insignificant thanks to the use of a narrow band- 
pass filter (chromatic refraction) . As the ADR is dependent 
on the zenithal distance, it may be significant even for a 
monochromati c light. So, the real reas on why ADR is negli- 
gible in case of lNeuhauser et all (|2006l ) is likely the geometry 
of their binary. Nevertheless, it is true that using wide-band 
filters makes ADR harder to calculate due to it's chromatic 
chara cter and, for example, stars' different colors (H elminiakl 
120091 ). 

Unfortunately, we did not collect any weather readings 
during our observing runs, so we had to use the canonical 
values of temperature and pressure and assumed 50 per cent 
humidity. For the Keck observations we assumed 2 times 
smaller pressure and the temperature of 0°C. This means 
that the real uncertainties of measured separations and po- 
sition angles are higher than the precisions given in Table 
[U To correc t for ADR, we us ed the semi-full approach as 
described by iHelminiakl (|2009l ). 

In order to estimate the maximum error due to ADR, 
we took the largest possible separation in our sample - 30.8 
as for GJ 873 1-2 (from Table©. We assumed that the dif- 
ference in the zenithal angles equals to the separation and 
the binary is seen 30 degrees above the horizon. In such 
an improbable case the maximum contribution to the er- 
ror budget coming from the temperature is 4 mas (if the 
real temperature were T = 230 K and pressure p = 1013.25 
hPa) and the fraction coming from the air pressure is smaller 
than 8 mas (for p = 613 hPa and T = 230 K). One should 
add that the bigger is the part coming from the tempera- 
ture, the smaller is the contribution from the pressure. So 
we may conclude that in this improbable case 8 mas is the 
maximum error and in most (if not all) our real observa- 
tions the uncertainty caused by the weather conditions is 
smaller than several mas. In the case of binaries observed 
with Keck II, the maximum error should be much smaller 
because the separations (and z differences) are smaller and 
other, more probable weather conditions were assumed. The 
maximum uncertainty scales linearly with air pressure and 
almost linearly with separation (maximum z difference) and 
temperature. In Section 5.1 we estimate in yet another way 
the observed rms of our astrometric measurements. 



4.3 Chip geometry and orientation 

Detectors are not perfectly rectangular, flat and perpendicu- 
lar to the light-path. At the astrometric precision necessary 



Table 4. Average pixel scales and North orientations for Hale 
telescope. 



Night 


Scale 


North 




(mas/pix) 


(degrees) 


23 Apr. (narrow) 


25.168(34) 


334.043(99) 


23 Apr. (wide) 


40.00(2) 


334.043(99) 


23 Jun. 


25.168(34) 


334.043(99) 


24 Jun. 


25.171(21) 


334.039(18) 


26 Jun. 


25.171(21) 


334.039(18) 


21 Aug. 


25.156(10) 


334.072(11) 


22 Aug. 


25.156(10) 


334.072(11) 


13 Nov. 


25.156(10) 


334.723(15) 



Table 5. The average separations between the stars 1 to 4 in the 
NGC 6871 open cluster. 



Pair 


p [mas] 


rms 


NGC 6871 1-2 


2171.719 


0.604 


1-3 


6906.003 


0.533 


1-4 


11303.461 


1.023 


2-3 


5257.245 


0.302 


2-4 


11463.147 


1.194 


3-4 


16144.625 


1.019 



to detect planets (well below 1 mas), one has to know how 
the camera's detector is distorted or how the pixel scale and 
the detector's orientation changes from epoch to epoch. For 
instruments mounted in the Cassegrain focus, as is the case 
for PHARO, the distortion changes with the telescope's po- 
sition due to gravity. This effect is not present or negligible 
in the case of NIRC2 which is located in a Nasmyth plat- 
form. The astrometric calibration and distortion models are 
available for both cameras. The calibration for PHARO is 
more complicated and includes not only the geometry and 
orientation of the detector itself but also the influence of 
telescope's position and til t of the ch i p rela tively to the light 
path. This is described by iMetchevI (|2006f1. The distortion 
of the NIRC2 camera was investigated during it's pre-ship 
testing and is described bv lThompson et al.l (|200ll n. 

The calibration we carried out included deriving 
the average plate-scale of the chip and the position of 
the North direction with respect to the Cassegrain ring 
(CR). The nominal values are 25.10 mas/pix for the 
39.91 mas/p ix for the wide fiel d pixel scale 

' 200 ll ). 



narrow 
and 335' 



8 for CR (Hay ward ct al 



As it was 



shown by iMetchev fc Hillenbrand! (|2004) . the real values 



are different from the nominal one and usually change 
from epoch to epoch. As a base for our calibration, 
we adopted the meas u remen ts from 23 June 2002 by 
IMetchev fc Hillenbrand! (|2004l ) which are: 25.168 ± 0.034 
mas/pix and 334.043±0?099. We chose 4 stars in NGC 6871, 
marked in Fig. [T] as 1-4 which we believe to be members 
of the cluster (their positions with respect to the star No. 
5 changed in a similar way), and using their relative po- 
sitions we have recalculated the average pixel scale in the 



3 See http://www.astro.ucla.edu/~metchev/ao.htmlj 

4 See http:/ /alamoana. keck. hawaii.edu/inst/nirc2/preship_testing.pdf 
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Figure 5. Separations of stars from the open cluster NGC 6871. The average value is plotted as a solid line and the corresponding rms 
in mas is shown. 



narrow field by assuming that their astrometric motion is 
not detectable. First of all, for every night and every of 
six possible pairs, we calculated a preliminary pixel scale 
and the North direction i ncorp orating the uncertainties from 
iMetchev fc Hillenbrand! (|2004j ). Later we averaged the re- 
sults for every single night. We checked if for two or more 
consecutive nights the pixel scale changed, and if it did not, 
we averaged the result over the number of nights. This pro- 
cedure allowed us to improve our pixel scale's uncertain- 
ties w ith respect to those given bv lMetchev fc Hillenbrand] 
l|2004 ) who used only one pair while we used up to 18 (six 
pairs, observed during 3 nights). An adequate procedure was 
carried out for the CR orientation angle. 

For 3 nights of August and November 2002 we obtained 
the value of 25.156 ± 0.010 mas/pix which is in poor agree- 
ment with the previous value but then the difference in 
separations of the stars was clearly seen. We also recalcu- 
lated new pixel scale for June 24 and obtained 25.171 ± 
0.021 mas/pix which shows that there was actually no scale 
change during one night. We also assumed the same plate- 
scale for the night of June 26. The cluster NGC 6871 was 
not observed in April 2002 but the results for GJ 458 sug- 
gest that the pixel scale and orientation was the same as 
in June. We assumed that the orbital period of GJ 458 is 
long enough that no orbital motion can be seen after 62 
days and the separation remains constant. Comparing this 
system with GJ 195 which is closer to the Sun and has a 
few times smaller angular separation, we may conclude that 
the peri od is much lo nger than 338 years (orbital period of 
GJ 195. iHeintdl 19741 ). probably close to 3200 years. We see 
only a small (however clear) motion in GJ 195 and we should 



expect at least a 10 times smaller movement in GJ 458 which 
would be below the detection. Thus for April 2002 we used 
the same pixel scale as for June 23 with its relatively high 
uncertainty. 

This binary was also observed in the wide field mode 
on 23 April at the beginning of the night. By combining the 
measurements from the wide and narrow fields, we obtained 
40.00±0.02 mas/pix as a pixel scale for the wide field and no- 
ticed no change in the CR angle. Between June and Novem- 
ber we noticed two changes of the CR orientation. The values 
of CR for August and November were: 334.072 ± 0.011 and 
334.723 ± 0.015 degrees respectively. The recalculated value 
for June 24 and 26 was 334.039 ± 0.018 degrees. We should 
also note that the position angles were computed from Ax 
and Ay counted in pixels along the chip's axes, so the bigger 
separation, the smaller uncertainty in 9. 

Values of the derived plate-scale and North positions 
are summarized in Table [4] For the stars in the open clus- 
ter NGC 6871, we show the average values of separations 
together with their rms' (Tab. [S] Fig. [SJ). The average val- 
ues are calculated from 5 epochs corrected for the new pixel 
scales. Such a calibration is not perfect and these plate-scale 
and North direction values are subject to possible system- 
atic errors e.g. due to a limited knowledge of the weather 
conditions. 

Examples of the distortion and the results of employing 
its model can be seen in Figure [6] We show the uncorrected 
and corrected measurements of separation (Ax, Ay) for two 
stars in the NGC 6871 field. The distortion is clearly seen 
in the x-axis (the top- left panel). It is worth noting that for 
PHARO the distortion is much bigger in x than in y, but 
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in the j/-axis the random scatter is about 50% bigger. The 
histograms demonstrate that after the correction, we are 
able to obtain a gaussian statistics (the middle panels) and 
the Allan variance (AV) shows no obvious signs of systematic 
errors (the bottom panels). 

The real average plate-scales of NIRC2 were 
found to be in agreement with the nominal values 
(IMetchev fc Hillenbrand! l2004h but the ?/-axis was r otated 
by 1?24 clockwise from North l|Metchev et al.ll2005T ). Un- 
fortunately, due to a small number of useful images in our 
data set, we were not able to perform proper tests and our 
own calibrations. Hopefully, the location of the camera on 
the Nasmyth platform grants its stability. In particular, the 
results for 56 Per and GJ 569 B where three different field 
rotator positions were used demonstrates that the precision 
of the field rotator of NIRC2 is better than 0? 1 during one 
night. 



5 ASTROMETRY 

The astrometric measurements are presented in Table [5] 
where for each pair of stars and epoch (MJD) the separation 
p \mas\ and the position angle 6 [deg] is given. The lcr errors 
were calculated using the gaussian statistic of Aa; and Ay. 
The plate-scale uncertainty is included into the separation 
error, but not into A9 for Hale observations. It is because 
we wanted to show how small changes can be noticed be- 
tween two nights where the same chip orientation is present 
(see GJ 860 1-2 in August - MJD=52509 and 52510). CR 
orientation uncertainties are about one order of magnitude 
bigger, so they would dominate the 9 error budget. The un- 
certainty in 9 for GJ 300 is also underestimated because this 
system was observed with only one position of the field ro- 
tator of NIRC2. This is however not the case for 56 Per and 
GJ 569 B where different values of the field rotator were 
used for dithering. Hence the resulting formal error of the 
position angle presumably more realistically reflects the ac- 
curacy of this AO system. For systems observed more than 
once, the orbital, parallactic and proper motion can be seen. 
Even for the long-period binary GJ 195 (P ~ 338 w; lHeint3 
1974) there is a clear signature of the orbital motion. Also 
a closer inspection of GJ 873 reveals a motion of the double 
secondary system (see Fig. [8}. GJ 873 B is probably a real 
binary but at a different distance from the Sun than GJ 873 
A (parallactic motion is present). 

For 5 of our binaries (GJ 195, GJ 352, GJ 569 B, GJ 661 
and GJ 860), the orbital solutions are known and can be 
found in the Washington Double Star Catalog (WDSfl The 
corresponding orbital elements are presented in Table[7] The 
quality of the orbit is represented by the parameter q - the 
smaller value, the more accurate the orbital solution. In all 
five cases the orbital solutions are not perfect but only for 
GJ 195 the elements are poor due to a long period of the 
binary. 

In Figure [7] the comparison between the orbit and our 
measurements is shown. There are several possible sources 
of the discrepancy between our data and the orbits: (1) the 

5 Sixth Catalog of Orbits of Visual Binary Stars 
http: / / ad. usno.navy.mil/wds / ort>6 / orb6frames.html 



quality of the orbits; (2) the uncertainty in the pixel scale; 
(3) an imperfect ADR correction. The level of the discrep- 
ancies (~ 100 mas for GJ 195 and GJ 860, ~ 20 mas for GJ 
352) seems to favor the first explanation. Also, if the cause 
was the ADR correction, one would expect a much higher 
scatter (see next section). 



5.1 Overnight and long term astrometric precision 

During one night, for most of the binaries observed with the 
Hale telescope we were able to go down below 500 pas in 
astrometric precision of p and in some cases below 100. As 
expected, the precision is better for objects for which more 
single images were obtained. For pairs with similar bright- 
ness of the components, the astromertic error is smaller than 
for pairs with a high brightness difference. This is due to the 
poor S/N of the faint component as well as a need not to 
saturate the bright one. It is imaginable that the usage of 
weighting might improve the precision a little. 

In some cases when the stars are located on two different 
parts of the chip's mosaic, the astrometric errors are larger. 
In particular, for a very frequently observed, relatively close 
binary GJ 661 we achieved the highest overnight precision of 
38 pas. Despite the fact that in this case the distortion cor- 
rection is not perfect due to a location of the binary around 
the centre of th e chip where all 4 parts of the mosaic meet. 
Metchcv (2006) suggest not to use this area because of small 
differences in the chip's components geometry. It seems pos- 
sible that for a binary like GJ 661, one may be able to achieve 
a precision even below our 38 pas in one night. 

We achieved a similar level of overnight precision with 
the Keck II/NIRC2. The main difference between the two 
data sets is the significantly lower number of images. It is 
quite surprising that with 10 frames taken with the 'narrow' 
camera we reached ~ 120 pas precision for GJ 300. Never- 
theless, one should treat this value cautiously. A low number 
of useful images does not allow for a particularly accurate 
calibration. 

Three systems GJ 661, GJ 860, GJ 873 and the open 
cluster NGC 6871 were observed more frequently than the 
remaining targets. Using their measurements, we can esti- 
mate the astrometric accuracy of the Hale telescope over a 
120-140 day time span. For the open cluster we take only 
pairs with star No. 5 which we believe is not a member of 
the cluster. Having 5 or 6 (for GJ 661) epochs, we can fit a 
2-nd order polynomial to the measured separations. Such a 
polynomial is sufficient to model the proper, paralactic and 
orbital motion of a close pair of stars. The fits are shown in 
Figure [8] The rms' for all fits are collected in Table [8] Note 
that again the resulting rms' are worse for pairs with very 
high brightness ratios such as those for GJ 860 with the star 
no. 3 and GJ 873 with the star no. 1. These rms' are also 
obviosuly higher than single night precisions for the corre- 
sponding pairs of stars and can be treated as en estimate for 
the true astrometric errors incorporating systematic effects 
due to an imperfect plate scale calibration, limited knowl- 
edge of the weather conditions and long term astrometric 
stability of the telescopes/cameras which could not be ac- 
counted for with our limited calibrations. 
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Table 6. Separations and position angles of investigated stars. 
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Table 7. Orbital elements of five of our binaries for which orbital solutions arc available from WDS. 



Star 


p [yr] 


a [mas] 


e 


i [°] 


q°] 


«["] 


t [MJD] 


WDS ID 


q 


GJ 195 


338. 


3720 


0.0 


65. 


168.5 


0.0 


55197 


05167+4600 


5 


GJ 352 


18.4 


630 


0.29 


143. 


48. 


285. 


45663 


09313-1329 


3 


GJ 569B 


2.424 


90.4 


0.312 


32.4 


321.3 


256.7 


51821 


14545+1606 


2 


GJ 661 
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image number image number 





Figure 6. An example of the distortion and its correction for two stars in NGC 6871 cluster. Left panels refer to X (a) component, right 
ones to Y (5). Top panels: the measurements before (+) and after (x, shifted) the distortion correction. Middle panels: the histograms of 
the measurements (bin width of 0.05 pix) with a gaussian fitted to the corrected measurements (left). Bottom panels: the Allan variance 
of the uncorrected (dot-dashed) and corrected (solid) measurements and an infinitely long, white-noise signal with a given in the middle 
panel (dotted line). 



5.2 Detection limits 

The astrometric signal, Q [fJ,as], of a planet with a semi- 
major axis a [AU] and mass Mp (Jupiter masses) in a cir- 
cular orbit around a star at a distance d (parsecs) and mass 
Ms (solar units) is given by (jPravdo fc Shaklanlll996l ) 



Obviously the same relation can be used for an S-type 
planel|f| in a wide binary system. Assuming that an astro- 
metric signal above 3<r p can be treated as a real one, B in 
the equation [4] may be replaced by 3cr p [mas] (from Table 
After changing d [pc] to parallax 7r [mas] we obtain 



(4) 



6 S-type, or satellite-type planet in binary /multip le is a planet 
orbiting only one of the components jDvoraklll984l) . 
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Figure 7. Relative orbits from WDS compared with our astrometry. Right pannels are zoomed in on the secondary's position. The 
formal error bars are smaller than the symbols used. The discrepancy is probably dominated by the quality of the orbital solutions. 
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Table 8. The rms of 2-nd order polynomial fits to the measure- 
ments of p for the most frequently observed objects 



Pair 


rms 


No. of 


time span 




[mas] 


nights 


[d] 


GJ 661 1-2 


0.282 


6 


122 


GJ 860 1-2 


0.216 


5 


143 


1-3 


2.154 






2-3 


2.595 






GJ 873 1-2 


1.127 


5 


113 


1-3 


1.282 






2-3 


0.309 






NGC 6871 1-5 


0.872 


5 


143 


2-5 


0.717 






3-5 


0.763 






1-5 


2.671 







aM P [AU ■ Mj] = 1562.5 



o p Ms 



In the above aMp = 4 [AU • Mj] means that we can detect 

1 Mj (or more massive) planet in 4 AU (or wider) orbit, or 

2 Mj planet in 2 AU orbit, etc. 



(5) 



The detection limits for the binaries, in which at least 
one stellar mass is known or can be estimated are collected 
in Table [9] Subscripts I and refer to the order of stars 
given in the first column. As one can see, in principle it 
is possible to achieve a sufficient astrometric precision to 
detect a massive planet or a brown dwarf with the Hale and 
Keck telescopes. As we have demonstrated in the previous 
section, currently a longer term precision is up to several 
times lower than the one achieved over one night. However, 
one can use the rms' from Table[8]with the Equation [5] and 
compute long-term planetary detection limits. These long- 
term limits' are also listed in Table [9] 
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Figure 8. A 2-nd order polynomial fits to the separations of stars from Table [8] 



6 CONCLUSIONS 

Nearby binary and multiple star systems are excellent tar- 
gets for astrometric searches for extrasolar planets thanks 
to their proximity and the availability of natural reference 
stars necessary for relative astrometry. In our study of 12 
visual binaries/multiples and 1 open cluster with the Hale 
and Keck II telescopes and their adpative optics facilities, 
we have demonstrated that over one night one is able to 
obtain an astrometric precision reaching ~40 /ias. Such a 
precison is sufficient to detect Jupiter mass planets around 
components of binary nad multiple stars. However, in order 
to turn the precision into a long term accuracy required to 
detect planets, one must be able to account for ADR and 



the plate-scale changes. The ADR correction requires ac- 
curate weather readings and the plate-scale changes must 
be carefully calibrated. In our attempt to account for both, 
we were able to achieve a long term (over a 140 days time 
span) accuracy ranging from 0.2 to 2.7 miliarseconds that is 
several times larger then the corresponding overnight pre- 
cision, but still allowing for detection of massive planets or 
brawn dwarfs. Since we have had limited means to carry out 
the calibrations, it is quite possible that a higher long-term 
accuarcy can be reached with the existing AO facilities. 
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Table 9. aM P limits. 



Pair 


7T 








MJD 


ajMp j 


ajjMp j j Ref. 




[mas] 


[M ] 


[M ] 


[mas] 




[AU ■ Mj] 


[AU ■ Mj] 




GJ 195 1-2 


72.0 


0.53 


0.19 


0.31 


52509 


3.56 


1.28 1 










0.23 


52510 


2.65 


0.95 










0.12 


52592 


1.38 


0.50 


GJ 195 1-3 




0.53 




0.50 


52509 


8.28 


1 










0.44 


52510 


5.06 












0.24 


52592 


2.76 




GJ 195 1-4 




0.53 




0.54 


52509 


6.21 


1 










0.47 


52510 


5.41 












0.23 


52592 


2.64 




GJ 195 2-3 




0.19 




0.43 


52509 


1.77 


1 










0.46 


52510 


1.90 












0.24 


52592 


0.99 




GJ 195 2-4 




0.19 




0.51 


52509 


2.10 


1 










0.47 


52510 


1.94 












0.21 


52592 


0.87 




GJ 352 1-2 


94.95 


0.44 


0.41 


1.11 


52389 


8.04 


7.49 2 


GJ 458 1-2 


65.29 


0.40 


0.37 


0.40 


52389 


3.83 


3.54 3 










0.28 


52450 


2.68 


2.48 










0.36 


52451 


3.45 


3.19 


GJ 507 1-2 


75.96 


0.46 


0.37 


0.45 


52389 


3.70 


3.43 3 










0.33 


52450 


3.12 


2.51 










0.36 


52451 


6.25 


5.03 


GJ 569B 1-2 


101.91 


0.071 


0.054 


0.11 


52337 


0.012 


0.009 4 


GJ 661 1-2 


158.17 


0.379 


0.34 


0.079 


52389 


0.30 


0.29 5 










0.170 


52450 


0.63 


0.62 










0.044 


52451 


0.17 


0.16 










0.038 


52454 


0.16 


0.15 










0.058 


52509 


0.22 


0.21 










0.041 


52510 


0.17 


0.16 










0.282 


TOTAL 


1.05 


1.03 


GJ 767 1-2 


74.9 


0.44 


(J. 4 


0.17 


52509 


1.56 


1.42 3 










0.12 


52510 


1.10 


1.00 










0.09 


52592 


0.83 


0.75 


GJ 767 1-3 




0.44 




0.52 


52509 


4.77 


3 










0.28 


52510 


2.57 












0.39 


52592 


3.58 




GJ 767 2-3 




0.4 




0.43 


52509 


3.59 


3 










0.22 


52510 


1.83 












0.35 


52592 


2.92 




GJ 860 1-2 


249.53 


0.34 


0.2711 


0.078 


52450 


0.17 


0.13 5 










0.048 


52451 


0.10 


0.09 










0.061 


52509 


0.13 


0.10 










0.061 


52510 


0.13 


0.10 










0.087 


52592 


0.18 


0.15 










0.216 


TOTAL 


0.45 


0.37 


GJ 860 1-3 




0.34 




0.73 


52450 


1.55 


— 5 










0.41 


52451 


0.86 












0.94 


52509 


2.01 












1.37 


52510 


2.95 












1.34 


52592 


2.86 












2.154 


TOTAL 


4.60 




GJ 860 2-3 




0.2711 




0.75 


52450 


1.28 


5 










0.42 


52451 


0.71 












0.91 


52509 


1.29 












1.19 


52510 


2.02 












1.33 


52592 


2.54 












2.595 


TOTAL 


4.96 




GJ 873 1-2 




0.36 




0.57 


52450 


1.62 


3 










1.07 


52451 


3.04 












0.73 


52509 


2.07 












0.60 


52510 


1.71 












1.12 


52592 


3.20 












1.127 


TOTAL 


3.29 




GJ 873 1-3 


198.07 


0.36 




0.58 


52450 


1.65 


3 










0.85 


52451 


2.41 












0.76 


52509 


2.16 












0.61 


52510 


1.73 












1.03 


52592 


2.93 












1.282 


TOTAL 


3.65 




GJ 9071 1-2 


72 


0.53 


0.49 


0.20 


52509 


2.22 


2.05 3 










0.26 


52510 


2.89 


2.67 










0.24 


52592 


2.66 


2.46 



References: (1) Fisher & Marcy 1992; (2) Sodcrhjclm 1999; (3) Harmanec 1988; (4) Zapatero Osorio et. al 2004; (5) Delfosse et al. 2000. 

Note: If MJD is 'TOTAL', the limit refers to the rms of the fit given in Table — an estimate of a long-term astromctric precision for a given pair of stars. [8~1 
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Figure 8 - continued 
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